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Abstract: We extend a recent scenario of Kachru, Kallosh, Linde and Trivedi to fix the 
string moduli fields by using a combination of fluxes and non-perturbative superpotentials, 
leading to de Sitter vacua. In our scenario the non-perturbative superpotential is taken 
to be either, the racetrack scenario or the M = 1* superpotential for an SU(A^) theory, 
originally computed by Dorey and recently rederived using the techniques of Dijkgraaf-Vafa. ^ — ^ 
The fact that this superpotential includes the full instanton contribution gives rise to the 
existence of a large number of minima, increasing with N. In the absence of supersymmetry 
breaking these correspond to supersymmetric anti de Sitter vacua. The introduction of 
antibranes lifts the minima to a chain of (non-supersymmetric) de Sitter minima with 
the value of the cosmological constant decreasing with increasing compactification scale. 
Surprisingly a similar picture occurs for the simpler system of the racetrack scenario. The 
relative semiclassical stability of these vacua is studied. Possible cosmological implications 
of these potentials are also discussed. 
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1. Introduction 

Supersymmetry breaking and moduli fixing have been the main obstacles for string theory 
to make contact with low energy physics. These questions are also essential for the study 
of any possible cosmological implications of the theory. Over the years there have been 
several proposals to solve these problems. Supersymmetry breaking mechanisms include 
the effect of non-perturbative field theoretical effects such as gaugino condensation and, 
more recently, the explicit breaking due to the presence of antibranes or brane intersections 



in low scale string models. The remnant potential for the moduli fields in these cases is 
such that the global minimum is either anti de Sitter space with a very large (negative) 
vacuum energy or the potential is of the runaway type, towards infinite extra dimension 
and/or zero string coupling. The problem is then not how to break supersymmetry but 
actually what to do, after breaking it, with the remaining potential for the moduli. 

Recently there has been interesting progress towards the solution of the continuum 
vacuum degeneracy problem in string theory. The introduction of fiuxes of RR or NS-NS 
forms |l|~||5| has proven to be very efficient to fix many of the moduli fields, including the 
dilaton. However in typical orientifold (or F-theory) compactifications, the overall Kahler 
modulus cannot be fixed by this effect [^, ^. 

An interesting proposal to also fix this modulus, due to Kachru, Kallosh, Linde and 
Trivedi (KKLT) Q, combines the fluxes with non-perturbative superpotentials that have 
been discussed in string theory in the past, consisting of a single exponential of the cor- 
responding modulus. The modulus is then stabilised for a supersymmetric anti de Sitter 
point. The further inclusion of anti D3 branes breaks supersymmetry explicitly and can lift 
the minimum of the potential to a de Sitter minimum with varying value of the cosmological 
constant, depending on the different parameters of the theory. Even though this is a very 
simple set-up, involving much fine tuning and the simplest non-perturbative superpoten- 
tial, it represents a concrete step forward towards fixing the moduli after supersymmetry 
breaking, with a potentially realistic value of the vacuum energy. 

In a different direction, progress has also been made in the understanding of non- 
perturbative effects in supersymmetric field theories. In particular, several techniques have 
been developed that allow the computation of the exact non-perturbative superpotentials. 
In some simple cases they have been derived in closed form encoding the infinite sum of 
exponentials of the inverse gauge coupling expected from the full instanton and fractional 
instanton effects. 

In this note, we slightly extend the KKLT proposal by considering more general non- 
perturbative superpotentials than the single exponential considered in KKLT. In particular 
we consider the non-perturbative superpotential for a supersymmetric theory for which the 
exact superpotential has been computed explicitly, including the infinite instanton sum. 
This is the so-called M = 1* theory 0, §. 

The structure of the potentials generated in this way is such that there are many 
AdS supersymmetric minima in the absence of the antibranes and many de Sitter min- 
ima in their presence, with intermediate configurations having both dS and AdS vacua, 
all non-supersymmetric. The minima are such that the value of the vacuum energy de- 
creases with increasing value of the compactification scale. This rich vacuum structure 
may have interesting physical implications and can serve as a concrete nontrivial example 
illustrating the possible 'landscape' of string theory |^. It is also similar to the staircase 
potentials proposed by Abbott ||l^ and to the models presented in 12, 13 1 exhibit- 
ing a dynamical relaxation of the cosmological constant. Furthermore, there is enough 
freedom to fine tune one of the minima to have a cosmological constant as small as re- 
quired. The presence of several de Sitter minima could also lead to interesting realisations 
of infiation. 



We also consider simpler cases in which the superpotential is a finite sum of exponen- 



tials as it appears in the much studied racetrack scenarios [14|. These are the simplest 
extensions of the mechanism of KKLT and, in the cases that lead to several minima, they 
serve as simple examples to study the transition between different vacua. 

This article is organised as follows: After briefly reviewing the effect of fluxes to fix the 
moduli in the next section, we start considering in section ^, the general scalar potential for 
the Kahler modulus fleld for arbitrary superpotential. In section Q we flrst consider the case 
of superpotentials with two exponentials which can give rise to one or two different minima. 
Section ^ discusses the potential for the Af = 1* theory with its rich vacuum structure with 
and without breaking supersymmetry. Section |6| is dedicated to the stability analysis of 
the different minima. We recover in particular the results of KKLT about the life time of 
the dS minimum with smallest positive value of the cosmological constant which is much 
larger than the age of the universe but smaller than the Poincare recurrence time. We also 
discuss the relative probability for tunneling between different minima. 

2. Fluxes and moduli fixing 

Type-IIB strings have RR and NS-NS antisymmetric 3-form field strengths, and F3 
respectively, that can have a (quantised) flux on 3-cycles of the compactification manifold. 



42// F^ = M.. 



1 



47r2 



a Jb 



Hi = -K, (2.1) 



where K and M arbitrary integers and A and B different 3-cycles of the Calabi-Yau 
manifold. 

The inclusion of fluxes of RR and/or NS-NS forms in the compact space allows for the 
existence of warped metrics that can be computed in regions close to a conifold singularity 
of the Calabi-Yau manifold, with a warp factor that is exponentially suppressed, depending 
on the fluxes, as: 

Therefore fluxes can naturally generate a large hierarchy. Here Qs is the string coupling 
constant. 

Fluxes have also proven very efficient for fixing many of the string moduli, including 
the axion-dilaton field of type-IIB theory S = e'^ + ia. A very general analysis of orientifold 
models of type IIB, or its equivalent realisation in terms of F-theory, has been done by 
Kachru and collaborators [^, In the F theory approach, the geometrical picture corre- 
sponds to an elliptically fibered four-fold Calabi-Yau space Z with base space M. and the 
elliptic fiber corresponding to the axion-dilaton field S. 

The consistency condition in terms of tadpole cancellation implies a relationship be- 
tween the charges of D-branes, 0-planes and fluxes that can be written as follows: 

A^D3-iV53 + iVfl,. = ^, (2.3) 



where the left hand side counts the number of D3 branes and antibranes as well as the flux 
contribution to the RR charge: 



iVfiux = :^^ / H^AFs. (2.4) 

The r.h.s. of ( |2.3D refers to the Euler number of the four-fold manifold Z or in terms of 
orientifolds of type IIB, to the contribution of the D3-brane charge due to orientifold planes 
and D7-branes. Here kiq refers to the string scale in lOD and T3 to the tension of the D3 
branes. 

The fluxes generate a superpotential on the low-energy four-dimensional effective action 
of the Gukov-Vafa-Witten form 0: 



W= / GsAn, (2.5) 
Jm 

where G3 = — iSH-^ with S the dilaton field and is the unique (3, 0) form of the 
corresponding Calabi-Yau space. 

In the simplest models there will be one single Kahler structure modulus defining the 
overall size of the Calabi-Yau space which we denote by T = + i6 where r is the scale 
of the extra dimensions and b an axion field coming from the RR 4- form (T = ip in the 
conventions of [Q, §). The relevance of this modulus is that it is the one that cannot be 
fixed by the fluxes. Its Kahler potential is of the no-scale form, that is: 

K = K{ipi,ip*)-3\og{T + T*), (2.6) 

with K the Kahler potential for all the other fields <fi except for T. This implies that the 
supersymmetric scalar potential takes the form 

^SUSY = (K'~W.,WDjW) , (2.7) 



with K^^ the inverse of the Kahler metric K^j = didjK and DiW = diW + WdiK the 
Kahler covariant derivative. The T dependence of the Kahler potential is such that the 
contribution of T to the scalar potential cancels precisely the term — 3e^|M^p of the stan- 



dard supergravity potential, this is the special property of no-scale models ||T5[|. Since this 
potential is positive definite, the minimum lies at zero, with all the fields except for T fixed 
from the conditions D^W = 0. This minimum is supersymmetric if Dj-W = W = and 
not supersymmetric otherwise. 

Since the superpotential does not depend on T, we can see in this way that the fluxes 
can fix all moduli but T. In order to fix T KKLT procede as follows: 

1. Choose a vacuum in which supersymmetry would be broken by the T field, such that 
W = Wo^O. 



2. Consider a non-perturbative superpotential generated by euclidean D3-brane or by 
gaugino condensation due to a non-abelian sector of wrapped N D7-branes. For which 



the gauge couphng is = 27r^ = 27rReT. Which induces a superpotential of the 

9ym 

form Wnp = Be~'^'^'^/^ . Combining the two sources of superpotentials Wq + Wnp, 
they find an effective scalar potential with a non-trivial minimum at finite T and 
the standard runaway behaviour towards infinity, as usual. The non-trivial minimum 
corresponds to negative cosmological constant giving rise to a supersymmetric AdS 
vacuum. 

3. In order to obtain de Sitter vacua, KKLT, consider the effect of including anti D3 
branes, still satisfying the condition ( |2.3| ). This has the net effect of adding an extra 
(non-supersymmetric) term to the scalar potential of the form: 

with the constant D = 2aQT^/g^ paramterizing the lack of supersymmetry of the 
potential. Here ao is the warp factor at the location of the anti D3 branes and T3 the 
antibrane tension. The net effect of this is that for suitable values of D the original 
AdS minimum gets lifted to a dS one with broken supersymmetry. See figure 



Here we will modify the KKLT 
scenario in two ways. First, regard- 
ing the original fluxes, we can con- 
sider the supersymmetric configuration 
where Wq = with the form G3 be- 
ing of the (2, 1) type. That is we 
may include a non-vanishing Wo part 
in the superpotential but it will not 
be necessary. Second, regarding the 
non-perturbative superpotential, we ex- 
plore the simplest M = 1 supersym- 
metric model for which the full non- 




perturbative superpotential has been pjg^j,^ ^. scalar potential considered in § with 
computed, including the contribution one single de Sitter minimum, 
from all (infinite) instantons. This is 

the so-called M = 1* model, constructed from mass deformation of = 4 super Yang- 
Mills. 

Our first minor modification allows to start with an explicit supersymmetric model, 
before considering the low-energy non-perturbative effects. This avoids the need to fine 
tune the value of Wq in looking for non-trivial minima. Our second modification allows 
exploring the possibility of an exact instanton calculation, instead of a single instanton 
calculation as it is usually considered. We will see that this exact superpotential will have 
a constant piece (similar to Wq) and an infinite sum of exponential terms, allowing for a 
very rich vacuum structure. For completeness, we also considered simpler superpotentials 
including a sum of two exponentials as in the racetrack scenario. 



3. The general scalar potential 



3.1 The supersymmetric potential 

The standard J\f = 1 supergravity formula for the potential in Planck units reads 

VsvsY = e"^ K'WiWD~W - 3|VF|2 j , (3.1) 

where i,j runs over all moduli fields. As we already mentioned, K^^ = didjK where K 
is the corresponding Kahler potential and DiW = diW + {diK)W. In our case we are 
working with a model having only one Kahler modulus, (that is, h^'^{M) = 1) as we will 
be focusing on the T field. All other fields are assumed to have been fixed by the fluxes 
just as in |^] so the superpotential W will depend on the superfield T. 

Then our purpose is to study the scalar potential V{T) which also depends on the 
Kahler potential. We take the weak coupling result in 4-dimensional string models, namely 

K = -3log{T + T*) , (3.2) 

and neglect possible perturbative and non-perturbative corrections. For simplicity of no- 
tation we will write the field T in terms of its real and imaginary parts: 

T=X + iY. (3.3) 



Using (|3.1| ) and ( |3.2D the scalar potential turns out to be 

VsvsY = ^ {l\2XW' - 3W\^ - 3\W\^^ , (3.4) 

where by ' we understand derivatives with respect to the field T. To compute the super- 
symmetric minima of the scalar potential we need to calculate the derivative of (3.4) that 
is given by 

, _dV _ i2XW" - 2W'){2XW' - 3W)* - 2{W')*{2XW' - 3W) 
^ -fff- 2AX^ ' ^^-^^ 

Notice then that there can be two types of extrema. The supersymmetric extrema 
appear when 

2XW'-3W = {). (3.6) 

In this case the value of the potential at the extremum is clearly negative definite leading 
to an anti de Sitter vacuum. These extrema are minima provided that: 

\XW" -W'\>\W'\. (3.7) 

The nonsupersymmetric extrema occur at 

{XW" - W) = {W')*e^'^ , (3.8) 

where 7 = arg(2XW - 3W). 

In all cases considered here, our analysis shows that the condition (3^) is never fulfilled 
at the minima. Then all minima are supersymmetric and lead to negative cosmological 
constants. 



3.2 Supersymmetry breaking 

Following the lines of |^], to uplift the anti de Sitter vacua to de Sitter vacua we will 
introduce in our model a D3. This will break the supersymmetry of the system. The 
introduction of the antibrane does not introduce extra translational moduli as its position 
is fixed by the fluxes so it just contributes to the energy density of the system. This 



contribution is given by |16| 



5V = ^, (3.9) 

where the coefficient D is a function of the tension of the brane and of the warp factor 
oo and has the form D = where Qs denotes the string coupling. The coefficient 

D depends on the fluxes through the dependence of the warp factor ao on them and it is 
therefore quantised, although the range of values of the fluxes can be such that for practical 
purposes it may be considered as an almost continuum variation ^, 27|. 



If we add this supersymmetry breaking term to the scalar potential (^^) we have now 
that the expression for the scalar potential is given by 

V = VsvsY + 5V . (3.10) 

The effect of the supersymmetry breaking term depends on the range of values of the 
parameter D. If D is very small, the potential will not change substantially and the 
minima remain anti de Sitter. For a critical value of D one of the minima will move up to 
zero vacuum energy and then to de Sitter space. Continuing increasing D, more minima 
become de Sitter until all of them are de Sitter. If D is larger than another critical value 
the nonsupersymmetric term starts dominating the potential and starts eliminating the 
different extrema to make the potential runaway with X. The range of values of D will 
depend very much on the form of W and we will illustrate the behaviour in the examples 
below. 



4. Superpotentials with two exponentials 

Our main interest in this paper is to consider non trivial superpotentials including all 
instanton contributions in terms of an infinite sum of exponentials. For simplicity we 
will first discuss the simpler case of superpotentials with just two exponentials. We then 
consider the following superpotential 

W = Wo + Ae-"^ + Be-^'^ . (4.1) 

4.1 The racetrack scenario 

This superpotential has an interest by itself since it includes the standard racetrack scenario 
(when Wq = 0), very much discussed in order to fix the dilaton field at weak coupling. The 
origin of the so-called racetrack scenario is the condensation of gauginos of a product of 
gauge groups. For an SU(A^) x SU(M) group we would have a = 2it/N and h = 27r/M. For 
large values of A'^ and M, with M close to N the minimum is guaranteed to lie in the large 
field region. A similar result will happen for the field T taken to be the gauge coupling at a 



Figure 2: Contour plot for a racetrack type 
potential with several minima. 



Figure 3: Graph of the scalar potential with 
a = 0.1, A^1,B^ 5-104 and H/q = -lO^^. 



product gauge group coming from the D7 sector. Turning on Wq it is possible to find more 
than one minimum that, following the KKLT method, would then lead to many (order 10) 
de Sitter vacua depending on the values of the parameters Wo, N, M, A,B. In figure |2| we 
present a contour plot of the potential illustrating this behaviour with Wq ~ — 10~^ and 
all other constants in the range 1, 10. 

4.2 A finite instanton sum 

The potential takes a manageable form if we assume b = 2a which allows more than one 
minima and we will follow mostly for illustration purposes. In this case using (3.1), (^ 
and ( [4. ID the scalar potential turns out to be 

VsUSY 



6X2 



6B e^^^Wo cos{2aY) + 6B^ + 3A^ e^"^ + AaB^X + 
+ aA^X e^"^ + A e"^ (3 e'^^'^Wo + B{9 + AaX)) cos{aY) 



(4.2) 



Notice that this potential is periodic in Y with period In/a, and is also invariant under 
the reflexion Y — > —Y. 

From ( [4.2| ) it is easy to see that if Wq is negative (Wq < 0) the minima of the potential 
will be located at y = y^") = vrn/a, n = 0, ±1, ±2, . . . For each value of y*^") we will find 
a minimum at X = X^^\ Since the potential (O) is periodic in Y with period 27r/o we 
will just consider the cases n = and n = 1, being the rest of the cases copies of these two. 
We will denote these minima as {X^^\Y^^^), where + will denote the case n = and — 
will denote the case n = 1. In figure |3| we show an explicit example of the potential ( |4.2| ). 

At the minima of the potential we have that 

-2aX(±) - 

dxV\Y=Y(i)=0 =^ Wo = (i^e'*^* '(3 + 2aAW) + S(3 + 4aAW)j . 

(4.3) 



O 
O 

00 

O 
Cn 
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Using this we can get the values of the potential at the minima. Those values are given by 

= - ^ (±Ae-^'"' + 2i?e-^"^'"V . (4.4) 

-(+) ^ T/(-) 



Note that from (gj) we find that X^+^ > X(-) and then from (gj) V^'^' > '. Also 



note from (4.4) that we will always have two non-degenerate minima with null or negative 
value of the potential, that is, we will have either Minkowski or AdS vacua. 

On the other hand, if Wq is positive {Wo > 0), the points located at y = Y^"'^ = vrn/a, 
n = 0, ±1, ±2, . . . are no longer minima, but maxima. In this case we find that the scalar 
potential will have two degenerate minima in every In/a period of Y. Those minima fulfill 

_)_ y(2) _ 27r/a and X^^^ = X^'^^ (the reason for that is the invariance of (|4^ ) under 
Y — > —Y). This means that we should need more exponential terms in ( |4.2| ) in order 
to have at least two non-degenerate minima. Therefore we will concentrate in the case 
Wo < 0, as is the simplest case of that type. The analysis of the case Wq > is analogous 
to the case considered in Q, where they analyse the case of a scalar potential with just 
one minimum. 

Now, we will study what is the effect of breaking supersymmetry in this model. The 
introduction of the supersymmetry breaking term ( |3.9| ) acts in the following way: when 
D/[X^^^)^ <C V^^^'' the potential is almost not affected, and we still have two AdS minima. 
If we increase the value of D such that D/{X^^'^)'^ ~ Vq^\ the value of the scalar potential 
at the minima increases with D, so the minima will eventually become positive (dS minima). 
For larger values of D the minima become saddle points and then disappear. 

Also notice that the supersymmetry breaking term ( |3.9D does not involve Y, therefore 
the new potential will also be periodic in Y with the same period as before, and its minima 
will also be located atY = y(") = vrn/a, n = 0, ±1, ±2, ... A gain, for each value of we 
will find a minimum in X but now at different values than in the supesymmetric case. We 
will denote these values as X^"'\ For each of those values we will find again two different 
values of the potential. It is interesting to note that the position of the minima in the non- 
supersymmetric case {X = X(±),y = yW) does not vary significantly from the position 
of the minima in the supersymmetric case {X = X^^\ Y = Y^^^). 

Therefore, for a given value of D such that D/{X'^^^f ~ V^^^ the non-susy scalar 
potential will have two different minima that, depending on the value of the constants, can 
be either positive or negative. If we want to describe the present stage of the acceleration 
of the universe within this framework, we would like both minima to be dS minima, such 
that Vj+^ > ^0^"^ ~ 10-120 Planck units. This can always be done in this model for 
example by fine-tuning the value oi D. In figure ^ we show an explicit example where the 
two AdS minima that appear in the supersymmetric case shown in figure ^ become two dS 
minima in the non-super symmetric case. 

5. A/" = 1* theory 



In this section, we would like to study more general superpotentials were the full instanton 
sum has been computed, being the simplest such case the A/" = 1* theory. The J\f = 1* 



100 150 200 250 300 

(a) Cross section of the scalar potential a,tY = 
and Y = Tr/a. 

Figure 4: Example of a configuration with a 
r' = 3.5-10~^. In this case we find a configurat 
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(6) Cross section of the scalar potential along 
a line that joins the two minima. 

= 0.1, A = I, B ^ 5 ■ 10^ Wo = -10--\ and 
m with two de Sitter minima. 



theory is mass deformed = 4 super Yang-Mills in which the M = 1 chiral multiplets 
inside the = 4 Yang- Mills multiplet are given a mass. At the moment we do not 
have a concrete example where this theory appears at low energies on the D7-branes. We 
may think of ways that it could arise, for instance since the D-branes break one half of 
the supersymmetries, they tend to have a AA = 4 gauge theory inside. Having a model 
that breaks supersymmetry to M = 1 outside the brane would naturally induce masses 
to the chiral multiplets. However here we will take the M = 1* superpotential only as 
an illustrative example of what we can expect from full instanton contributions to the 
superpotential. 

5.1 Background 

Let us first review what the A^ = 1* theory is, starting from SU(iV), AA = 4 super Yang- 
Mills. The SU(A^), AA = 4 super Yang-Mills theory can be written in terms of AA = 1 
superfields as a gauge theory with three massless chiral superfields in the adjoint of 
SU(A^) with superpotential 

W = eijkTi<^i[<^j,<^k]- (5.1) 
A deformation of this theory by adding mass terms to these superfields 

AW = miTT^f, (5.2) 

with all mi ^ breaks supersymmetry to AA = 1. This is the M = 1* theory. There are 



further deformations of the original AA = 4 theory that may be considered [18|. 

The classical vacua of this theory can be found by solving OWx/d^i = for Wt = 
W + AW, which leads to 

[<^^,^j] = eijkmk^k (5.3) 

and therefore the solutions correspond to A^-dimensional representations of the SU(2) al- 
gebra, giving rise to the different phases of the theory. 

The massive (Higgs and confining) phases of this theory are well understood. They 
are labelled by a triplet of integers (p, q, k) with pq = N and < k < q. These phases 



interpolate from the full confining phase q = N to the full Higgs phase p = N. The 
exact superpotential for M = 1* was derived by using instanton techniques for the theory 
compactified on a circle. The compactification to three dimensions is a computational trick 
and it turns out that the superpotential is independent of the compactification radius Q . 
After integrating out the gauge fields, the superpotential depends only on the (complex) 
gauge coupling r and takes the form:^ 

%,,,fc(r) = E2{t) - + -\ . (5.4) 

Here and later in this section by £"„ we will denote the Eisenstein modular functions. 

The modular properties of the superpotential show that under a SL(2, Z) transforma- 
tion 

ar + b 

T (5.5) 

the effective theory in a given phase is not invariant but it exchanges different phases by 
changing the values of p, q, k. For instance r — > r + 1 requires (p, q, k) — > (p, q,k + p) and 
r — > — 1/t implies {p,q,k) (a,^,k') with a = gcd{q,k) and k' having a complicated 
dependence on p, q, k with k' = ii k = 0, indicating the exchange of Higgs and confining 
phases under this transformation. 

It is interesting to note that the validity of this exact non-perturbative result has 
been checked using string theory techniques [H| and more recently from Dijkgraaf-Vafa 



techniques 18], making this expression very robust. 



5.2 The scalar potential 



Following we will now promote the parameter r to a full M = 1 superfield that we 
identify with the modulus field^ T = —it. Then the superpotential ( |5.4D can be written in 
terms of T as 

Wp,g,k{T) = E2{T) - - i-\ , (5.6) 

q \q qJ 

The superpotential Wp^q^k{T) transforms as a modular form of weight 2 once the values 
of p, q, k are transformed accordingly. As we already mentioned, the net effect of a SL(2, Z) 
transformation is that it changes one massive phase to a different phase. 

For concreteness we will concentrate in the confining phase p = 1, q = N, for which 
it is enough to set A: = 0, since k ^ can be reached by a translation T ^ T — ik. In 



^An overall factor of order m'^ = ■^mim2m3 appears multiplying the superpotential, where rui are the 



24 

masses of the chiral superfields of the — 4 theory. This will scale the scalar potential by a particular 
scale, which we take as unity for simplicity, but need to keep it in mind when combining with the non- 
supersymmetric case and to take care of actual number estimates regarding the value of the cosmological 
constant. For consistency we have to take the scale set by the m^'s to be hierarchically smaller than the 
string scale. It is not yet clear if these small mass parameters will be achieved in explicit models. 

^In reference jl^] the corresponding field was the dilaton S rather than T then the difference between 
the potential calculated there and the supersymmetric potential we consider here is only in the factor of 3 
in the Kahler potential. The main difference appears when we consider the nonsupersymmetric case. 



this case the scalar potential will be given by (| 
derivative take the form 



W{T) 
W'{T) 



1 /T 

E2(T) E2 — 



U{E,{T)-E2{Tf) 



where now the superpotential and its 

(5.7) 

(5.8) 
2^^^ that 




We will mostly work with the En{T) expansions in terms of the variable q = e 
are given by 

00 00 
= l-24^ai(n)g", = 1 + 240 ^ ^3(71)9" , (5.9) 

n=l n=l 

where (7p(n) is the sum of the p*^^ powers of all divisors of n. E4 is a modular form of 
weight 4 so i?4(l/T) = T^E4{T), while E2 fails to be a modular form of weight two since 



E2i^^'^=-T'E2{T) + ^. (5.10) 

Since there is already a constant term in the expansion of the superpotential we will 
perform most of our analysis considering Wq = which means no supersymmetry breaking 
by the fluxes themselves. Contrary to the case with one single exponential, we will find 
many minima without needing to tune the value of Wq to obtain large compactification 
volume naturally once we break supersymmetry. We will see later how Wq affects our 
results. 

5.2.1 Numerical analysis 

In order to perform reliable computations with the Eisenstein series we use the 'weak 
coupling' or large radius expansion of W(T) in eq. ( |5.6| ) only when In X > N. For other 
ranges we can use the property ( 5.10| ). For instance, for 1 < 27rX < we can transform 
the E2(T/N) term to obtain 



W{T) = E2{T) + ^^E2(^)-:^. 



(5.11) 



Similarly, when 2-k X < 1^ we can trans- 
form both terms in W(T) to obtain 

T 

(5.12) 



W{T) 



1 N 
-j^E2{T) + ^E2 



N 


2 3 4 


X 


1.46 1.70 1.82 




-1.73 • 10~2 -1.60 • 10~^ -1.32 • 10^2 



Table 1: Minima of the scalar potential. 

For < 4 we have that the minima appear when 2tt X > N , so that one has to use 
the expression ( ^.7D for the superpotential. One finds the supersymmetric minima are at 
Y = N/2 and X (when < 4) at the values given in table ||. 

For > 5 the fiat direction at y = A^/2 turns into a saddle point and pairs of 
supersymmetric minima Tji, Ti2 on either side in the Y direction appear, such that Yn + 
Yi2 = N and Xn = Xi2. Examples of these are given in table |^. 

^We have to keep in mind that in these regimes there may be substantial corrections to the Kahler 
potential which are not under control. 
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6 


7 


8 
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X 


1.48 


1.65 


1.77 


1.85 


1.88 


Y 


2.24 


2.40 


2.63 


2.88 


3.13 




-1.35 ■ 10-2 


-1.37-10-2 


-1.27- 10-2 


-1.16 • 10-2 


-1.07- 10-2 



Table 2: Minima of the scalar potential for several cases. 
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Table 3: Minima of the scalar potential for N ~ 30, 40, 50. 





Figure 5: Graph for the scalar potential 
with N = 100. 



Figure 6: Contour graph for the scalar po- 
tential with = 100. 



For > 10 we have that the minima begin to appear when 2tt X < N, so that one has 
to use the expression ( 5.11 ) for the superpotential to compute the scalar potential (0). 
This needs to be done in order to have a convergent series expansion in (|5.9| ). Using this 
we find a lot of super symmetric minima in the scalar potential, noting that the number of 
minima increases with N. We list some of the results table ^. 

Also we present in figure |5| a 3D plot of the potential for the case N = 100 illustrating 
the rich structure of the potential with many supersymmetric AdS minima. In figure ^ we 
present a contour plot showing the location of the minima and how the landscape changes 
in field space. 



N 


Number of 
Minima 


Mine 

X Y Vmin 


Minx 

X Y Vtnin 


10 


2 


1.28 3.59 -1.11 • 10~^ 


1.85 3.32 -1.01 - 10-2 


20 


3 


1.76 7.84 -1.43 • 10"^ 


1.96 4.65 -6.92 - 10"^ 


30 


5 


1.57 8.40 -1.70-10"^ 


2.08 5.69 -5.21 - 10"^ 


40 


6 


1.44 8.73 -1.69 • 10"^ 


1.28 6.95 -5.73 - 10"^ 


50 


7 


1.59 10.98 -1.70 • 10"^ 


1.10 7.52 -3.46 - 10-^ 


60 


10 


1.45 10.78 -1.70 • 10~^ 


1.20 7.62 -1.62 - 10-^ 


70 


10 


1.53 12.62 -1.76-10"^ 


2.32 8.61 -2.67-10-^ 


80 


10 


1.62 14.45 -1.69 • 10"^ 


1.13 8.85 -1.07-10"^ 


90 


10 


1.49 13.75 -1.78 • 10"^ 


2.38 9.73 -2.16 - 10"^ 


100 


10 


1.55 15.28 -1.76-10"^ 


2.41 10.25 -1.97-10"^ 


200 


15 


1.52 20.98 -1.76 • 10~^ 


2.56 14.37 -1.06-10"^ 


300 


19 


1.60 26.00 -1.74-10"^ 


2.70 17.50 -7.26 - 10~* 


400 


20 


1.49 29.59 -1.80 • 10~^ 


2.68 20.20 -5.53 - lO"* 


500 


25 


1.47 32.20 —1.77 • 10 


2.70 22.60 —4.46 - 10 


600 


25 


1.50 36.31 -1.79 -10-2 


2.73 24.70 -3.74 - 10-^ 


700 


27 


1.47 37.80 -1.75-10-2 


2.75 26.60 -3.22 - 10-^ 


800 


31 


1.62 45.70 -1.72-10-2 


2.77 28.43 -2.83 - 10-^ 


900 


34 


1.56 46.10 -1.76-10-2 


2.78 30.10 -2.52 - 10-^ 


1000 


38 


1.48 46.50 -1.72 - 10-^ 


2.78 31.80 -2.28 - 10-^ 



Table 4: Number of minima for different values of N in the supersymmetric case. 

To illustrate the fact that we find an increasing number of minima when we increase N ^ 
and also other interesting facts, we present in table |^ some other results of the numerical 
analysis. Note that by Mine in table ^ we denote the global minimum of the potential and 
by Minx we denote the minimum with a largest (finite) value of X. 

From the information shown in table Q we may extract some general remarks regarding 
the minima. First, the number of minima grows quite fast when grows. Also the potential 
at the minimum with larger value of X increases with N . Furthermore, even though we 
do not have a general proof, for all the cases we have analysed, all the vacua of this theory 
happen to be supersymmetric. 

Finally, although it is not necessary in our case, following we have studied the 
effect of turning on Wq. Notice that with Wq = we have many minima but all of them 
correspond to not too large values of the compactification scale X. We have found that 
tuning Wq has the effect of allowing minima with very large values of X improving the 
validity of the field theoretical analysis of the potential. Otherwise the general structure 
of the potential remains the same. 

5.2.2 Analytical considerations 

The scalar potential (|3.4| ) appears to be, in general terms, too complicated to do an analyt- 
ical study of its minima for a superpotential such as the N = \* one. In fact, introducing 
the superpotential ( ^.11 ) into the condition ( |3.8D that is used to find the supersymmetric 
minima we arrive to non linear equations that cannot be solved analytically. 



In any case, some interesting features appear when we study the behaviour of the 
potential in the hmit {2tt X) N ^ + y^. In that hmit the superpotential W can be well 
approximated by the function 

N 6 

W r^l + —r-—. (5.13) 



This approximation is derived from ( 5.11 ) just by keeping the constant terms in the E2 
expansions and neglecting the exponential terms. 

Using the approximate superpotential in ( 5.13| ), we find that the condition for a super- 



symmetric minimum, i.e., 2XW' — 3W = 0, reduces to a cubic equation in X which can 
be solved analytically. This cubic equation is given by: 

^— + 2NT?j X - 168X2 + 2ATrX^ = 18N , (5.14) 

with Y given in terms of X by the following relation: 



y = ^A. + ^P=^. (5.15) 

The full form of the solutions for X and Y is not very enlightening, but they allow an 
expansion in powers of 1/A^ that happens to be more useful.^ 

9 3240 1 5015520 1 
Xmin = 3-^ H 5 — ^ H (5.16) 

r- / 99 1 502281 1 \ , , 

It is clear from these expressions that when N is large the solution for X will tend 
asymptotically to 9/7r, and Y will tend to ^/N. It is also possible to compute the value of 
the potential in this minimum (also as an expansion in power series of N) substituting the 
approximate expression for the superpotential W given in (|5T^ ) into dU). Its expansion 
in powers of 1/N is given by 

27r 1 50 1 , , 

^- = -27iV + 3^iV^ + ---- ^'-''^ 

In fact, this minimum is the minimum with the largest finite value of X. Also we have found 
that for N > 50 these resuhs (|5.14| ), ( |5.15D and ( |5.1^ ) agree well with those obtained in 



the numerical analysis (see table Q). This minimum is in fact the minimum that appears to 
have the largest value of the potential for large values of N. This will play an important role 
in the next section, where we explore the minima of the potential when the supersymmetry 
is broken. 



*The equation (5.14) has of course three solutions for X, but only one of the three solutions of the cubic 
equation gives a positive Xmin in the range of validity of the approximation. 



5.3 Supersymmetry breaking 



In this subsection we will study the changes produced in the structure of the vacua when 
we introduce in the scalar potential the supersymmetry breaking term (|3.S|). We will see 
that with the introduction of such a term it is possible to lift the vacua from anti de Sitter 
vacua to de Sitter vacua, for some range of the parameter D.^ 

5.3.1 Analytical considerations 

As in the previous case, the scalar potential ( p.l[l[ ) appears to be, in general terms, too 
complicated to do an analytical description of its minima. In fact, introducing the super- 
potential ( 5. lip into the scalar potential (3.10) and minimizing it to find the minima leads 



to non linear equations that cannot be solved analytically. 

In any case, similarly to the previous cases, some interesting features appear when 
one studies the behaviour of the potential in the limit 2t^ X N (-^)^ + O^Y'- As in 
the previous case, in that limit the superpotential W can be well approximated by the 
function (5.13). 



Using the approximate superpotential in ( 5.13| ), we find that the condition for minimum 



reduces to two polynomial equations m. X , Y . Those equations cannot be solved in general 
but for N large it is easy to find that Y is given in terms of X by the following relation: 



Using this relation ( 5.1S| ), we find a complicated polynomial equation on just X. This 
equation can be simplified for large to a cubic polynomial on X"^, and therefore can be 
solved analytically. In fact, what we found is that only two of the solutions are positive 
and then lead to two real positive solutions for X. Nevertheless just one of these two 
solutions is a minimum, the other one corresponding to a saddle point. The solution 
found is proportional to \/iV, with the proportionality constant given in terms of the 
supersymmetry breaking parameter D. That is, the solutions for the minimum are given 
hy X^^ = f{D)Ny^. 

The exact form of the function f{D) is a complicated expression, but in any case it is 
possible to extract some useful features from it. In fact, it is easy to see that f{D) is an 
increasing function on L>, until it reaches an upper bound Dm- For values D > Dm the 
values of the function f{D) are no longer real, this meaning that the potential would have 
no minima. Therefore, in the limit of validity of the analytical analysis, there is a bound 
for the supersymmetry breaking parameter D if we want the potential to have minima. 
This bound is given by 

^ < ^ (27 + 7^21) = Dm • (5.20) 
3dU 



^Recall that the A/" = 1* superpotential has a mass scale m^. This was irrelevant for the discussion of 
the supersymmetric case since it could be rescaled out of the potential, affecting only its absolute value at 
the minima. We have to keep this in mind when considering the range of variation of the parameters Wo 
and D which will carry now an extra scale determined by m. 



Therefore, for values of the supersymmetry breaking parameter D larger than this bound 
we will not find any minima of the scalar potential. 
When the bound in D is saturated we find that 



Y 12(-4 + V2I) ' ^ ^ 

These values for X and D allow an expansion in powers of 1/N of the potential in the 
minimum. We find that the value of the potential is positive and is given by 

= 0-0406^ + 0.0386^ + • • • (5.22) 

Also from the form of f{D) it is clear that one can fine-tune the supersymmetry 
breaking parameter D so that the value of the scalar potential in the minimum Xmin 
approaches arbitrarily to zero. In fact one can compute the value of D such that Vmin = 0. 
The full expression is not very useful, but it allows an expansion in power series in 1/A^, 
which is given by 

3 9 / 3 99 1 

20"20^V5iV + 1600^iV + '" ' ^^'^^^ 
and also the values of X^m and Ymm are 



^;i... = >/^iV + £-^\/^ + - (5.24) 



..0 /5 1 3463 / 1 

where the value of Y^^^ is obtained from ( p. 191 ). Also remember that V^jj^ = 0. 

In fact, we have found that for > 50 these results agree well with those obtained in 
the numerical analysis, as we will show in the following subsection. 

5.3.2 Numerical results 

In this subsection we present the results found in the numerical analysis of the scalar 
potential in the non-supersymmetric case. We have computed numerically the minima of 
the potential for D = Dq for several values of A^, such as A^ = 50, A^ = 100, A^ = 500 and 
A^ = 1000. The results obtained from the analysis are shown in table |5|. 

As a matter of comparison with the supersymmetric case, we show in figure |^ the 
scalar potential for A^ = 100 with supersymmetry broken. 

Finally we show in table ^ a similar information as the one shown in table Q in the last 
subsection, where we write the number of minima varying with A^. In table || we denote 
by Minx,Y the minima with largest value of X, Y. We can notice that the number of 
minima still increases with A^ but slower than in the supersymmetric case. An argument 
for this is that the non-supersymmetric term 5V tends to smooth out the potential as it 
moves the minima up, so some of the minima eventually disappear. Another implication of 
this term is that the minima tend to have larger value of the compactification scale, large 
enough to create a hierarchy as compared to the string scale, as desired phenomenologicaly. 
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Table 5: Minima of the scalar potential for several non-supersymmetric cases. 

Furthermore, unlike the supersymmetric case, the minima appear to be ordered in X, where 
by ordered we mean that the value of the potential decreases as X increases. The one with 
largest compactification scale has a smallest cosmological constant. 

In figure ^ we illustrate the effect 
of the non-supersymmetric term in the 
potential. The value of the potential 
at the minima is presented for several 
values of the parameter D. For D = Q 
we have the supersymmetric case with 
more minima and not ordered, the in- 
creasing of the parameter D will re- 
duce the number of minima, increase 
the value of the compactification scale 
at the minima and order them from 
larger to smaller values of the cosmo- 
logical constant as the radius increases. 

Turning on a nonvanishing value of Figure 7: Graph for the non-supersymmetric scalar 
Wo tends to change some of these con- potential with N = 100. 
elusions. The minimum with the largest 

value of X can be lifted to positive vacuum energy while some of the other ones remain 
anti de Sitter. Besides this, the number of minima stays of the same order. 

We can see that, in general, increasing the value of N substantially increases the 
number of vacua and the value of X at the minima. Notice that can be quite large 




AT 
iV 


NunibGr of 

\A in 1 in Q 

IViillilliCi 


Y Y 1/ ■ 


X Y V ■ 


in 


9 


O.t/O O.t/U u 


^ Q8 ^ Qn n 




9 


Q Q9 on 


Q Qo K nn n 


^n 


9 


4 49 1 


9 88 1 9 nn 7 fiQ . 1 n"~^ 




9 


4 qq 7 10 


■^40 ifi "^n 4fin.in~^ 




Q 
O 


\ 7 qzL n 


4 4fi 9n 48 '\'\(\ ■ ^ 


fin 

uu 


Q 
O 


^ RQ R 70 


4 nfi 91^ f^n no . in~4 


7n 




9Q Q 40 


4 17 "^n 1 n 9 49 . 1 n~4 


sn 

ou 




fi^ 1 no 


4 98 ^4 8n 1 Q9 . 1 n~4 


on 




7 04 1 70 


4 "^"^ 40 1 fiO . 10""^ 


1 nn 

-LUU 




7 Qp; 1 1 94 


4 R9 44 90 1 . lO"'^ 


9nn 


O 


1 09 1 00 


^ q4 qi ^^0 4 '^O . 10~^ 


■^nn 

ouu 


(J 


19 0*^ 1 q 40 


fi1 1 07 ^0 1 1 . 1 0~4 
U.Ul lUl.OU i.±o ±u 


4nn 


u 


1 Q 7A 99 40 


7 T'^ 1 94 ^0 7 ^9 . 1 0~^ 


5nn 


fi 
\j 


1 5 40 25 1 n n 


8 54 1 5fi 70 5 • 1 0^^ 


600 


8 


16.70 27.50 


9.29 211.00 4.21 • 10"^ 


700 


8 


17.97 29.70 


10.09 245.30 8.03 • 10"^ 


800 


9 


19.20 32.80 


10.81 253.80 2.79 • 10"^ 


900 


9 


20.30 33.70 


11.10 432.90 2.35 • 10"^ 


1000 


9 


21.34 35.51 


11.46 482.05 2.42 • 10"^ 



Table 6: Minima for different values of N in a non-supersymmetric case with D = Dq. 

in string theory and wc may attempt extracting results varying N over many orders of 
magnitude. However, there are limiations on the numerical analysis if N is too large, since 
many more terms in the series have to be considered and the accuracy of the results would 
not be guranteed. Also computer and woman time to perform the analysis would increase 
substantially. Notice also that the number of minima we present here is actTially a lower 
bound since, again, the accuracy of the numerical calculations could have missed some of 
them. In particular, minima which would be very close to each other could be missed. 

6. Stability of the vacua 

Having a situation with many nonsuper symmetric vacua with positive cosmological con- 
stant implies that the relative stability of this vacua should be considered. In particular if 
our own universe may be described by one of them, it is important to know the lifetime of 
each of the vacua against decay to any of the ones with smaller value of the cosmological 
constant and compare it with the age of the universe. 

In the general case with many vacua the stability question is very much model and 
vacuum dependent. Given any particular minimum of the potential there is a nonvanishing 
probability towards tunneling to any of the other minima. The lifetime will depend on the 
height of the potential and the side of the local maxima or saddle points that separate 









\MD = 


0.06 


\±iD = 


0.08 




0.1 


\MD = 


0.12 


[A]D = 


Do(= 0.139) 



Figure 8: Potential versus X for N ~ 100 and different values of D. 



two different minima. The general analysis is therefore very complicated. We will limit 
here to the well controlled case of two notrivial minima that we found in the case with two 
exponentials. This may serve as an illustration of the techniques used in the general case. 

We know that the effective scalar potential has the standard runaway behaviour to- 
wards infinity in the radial modulus X. The existence of the runaway vacuum at oo in field 
space with zero energy is a common feature of all string theories This means that 

none of the two minima that appear in the scalar potential (see figures |3|, will be stable, 
as they are unstable, not only to decay from one minimum to the other minimum, but also 
to decay to the runaway minimum at X ^ oo, by means of either quantum tunneling or 



thermal excitation over a barrier [21 



As we already mentioned, if we want to describe the present stage of the acceleration 
of the universe using a model of the kind presented here, we would like the last de Sitter 
minimum to be so that V ~ 10^^^'' in Planck units. This can always be done in these models 
just by adjusting the value of the different parameters of the theory. In the case presented 
in section ^ we can achieve this by fine-tuning Wq and D, while in the Af = 1* such a value 
for the cosmological constant can be obtained by slightly modifying the value of -D = Dq 
(for which the value of the potential is zero in the last minimum) by a small amount. 



Doing so one can always get the desired value for V in the last de Sitter minimum. We will 
restrict here to the analysis of the case considered in section |3|, so we will have two minima 
with different values of the scalar potential Vq^\ Vq ■* such that Vq^^ > ^ ~ 10"-'^^''. 
Therefore we will have a model with several decay possibilities. In fact from the minimum 
at Vq^^ it is possible to decay either to the minimum at V^'' ^ or to the minimum at X — > oo, 
while from the minimum at Vq ■* it is just possible to decay to the minimum at X — > oo. 
We will also comment some features of the M = 1* case. In order to analyse this issue we 
will review several features of tunneling theories taking into account gravitational effects, 



following the original work of Coleman and De Luccia [21|. 
6.1 Vacuum decay 

Let us consider a theory of a scalar field ip with a potential V{(p) which has two local 
minima at (fi, with V{'^i) > V{(p2)- Both of the minima are stable classically but the 
vacuum state at ip = ipi (that is, the false vacuum) is unstable against quantum tunneling 
and will finally decay into the true vacuum state at ip = ip2, this vacuum decay proceeding 
through the materialisation of a bubble of true vacuum within the false vacuum phase. 
The tunneling action is given by 

S{ip) = J d'x^g (^-\r + \{d^f + V{^)^ , (6.1) 

with a tunneling probability between two vacua, per unit time and unit volume, given by 

P{^) ^ = e-^(^)+S(^i) . (6.2) 

Here is a solution of the equations of motion, which is usually referred to as "the bounce" , 
and S{'pi) is the euclidean action in the initial configuration ip = (pi. In the limit of small 
energy density between the two vacua Coleman and De Luccia showed that the coefficient 
B can be calculated in a closed form. Although in the general case it is usually very 
difficult to find analytical solutions for the Coleman-De Luccia instantons and calculate 
the probability of decay through quantum tunneling, the computation can be simplified 
within the range of validity of the thin-wall approximation. This approximation is valid 
when the thickness of the transition region between the true and the false vacuum is small 
compared with the radius of the bubble. In Minkowski space, the condition T^nin <Si Vma.^ 
usually means that the thin-wall approximation is applicable, so the false vacuum state 
will decay through the materialisation of a bubble of true vacuum within the false vacuum 
phase, which is a quantum tunneling effect (by T^ax we denote the high of the de Sitter 
maximum that separates two minima). That condition is usually fulfilled in the model 
considered here, and we will explicitly check later that in those cases we are always within 
the limits of the thin-wall approximation. 

Now let us consider the case in which we have a vacuum with positive cosmological 
constant (a dS vacuum) and a vacuum with null cosmological constant (a Minkowski vac- 
uum).^ For those cases Coleman and De Luccia found that the decay will be produced 

''This is always our case unless in the decay from Vq'^' to Vq~\ that are both de Sitter vacua. Never- 
theless, as > yj"^ ~ 10-12°, it is a good approximation to consider that the decay is also from dS to 
Minkowski. 



through a nucleation of a bubble of radius^ p = ^^^^^^^ , where e denotes the energy density- 
difference between the two vacua e = V{(pi) — V{(p2)- As the final vacuum has null cos- 
mological constant V{ip2) = 0, then e = V{ipi). Also Si denotes the tension of the bubble 
wall, that is given by 

Si= r dip V^Viip) , (6.3) 
The coefficient B is then given by 

The thin- wall approximation is justified in this context if p and \/3/e are large compared 
with the range of variation of the scalar field (f, that is, p, \/3/e ^ Aip. 

6.2 Comparing tunneling probabilities 

In this subsection we will compute the probabilities of all possible tunneling trajectories 
between the different vacua of the potential ( |3.10| ), (|4.2| ). Nevertheless before doing so we 
need to make a remark: In our case we have a potential that depends on two scalar fields. 
When two or more fields are involved obtaining the bounce (p, and therefore computing 
tunneling probabilities, becomes a much more difficult task (see for example |2^), except 
in trivial cases. 

An example of a trivial case is to consider the tunneling from any of the minima 
to the runaway minimum at X ^ co. In those cases the lines Y = 0, vr/a are always 
minima of the potential in the Y direction, so we can consider a one-dimensional potential 
y(±,oo) _ Y(^x, Y = 0, vr/a) and analyse the tunneling probability in the standard way (we 
would obtain one-dimensional graphs like the ones shown in figure |^) . Also note that the 
scalar field (p is defined in such a way that the kinetic term in the action (^]^) is canonical. 
In our analysis, the complex scalar field T has a kinetic term that is given by the derivatives 
of the Kahler potential (|3^), and is given by ^x)^ {dXdX + dYdY). If we consider that 
Y is fixed then the kinetic term coming from the Kahler potential reads |(9 In X)^, so the 

canonical field would be of the form (^(^'°°) = In X. 

On the other hand, the case of computing the tunneling probability from the minimum 
^0^+^ to the minimum Vq is more subtle, as we cannot trivially reduce it to the one- 
dimensional case in the same way as the previous case. It is however possible to find 
a lower bound for the tunneling probability by replacing the multi-field potential by a 
suitably chosen single field potential, as any chosen trajectory will have larger tension ( |6.3| ) 
than the minimum one. As we are interested in comparing probabilities such a bound 
will be enough for us. Therefore we can consider that the line that joins the two minima 
together is a good approximation for the bounce (as we already mentioned, this is an 
upper bound for the tension, as the bounce is the one that minimises the action). Then 
we can consider that Y = ^^(+'.-) ^ + ^ where by AX, we denote the difference 

in the value of X,Y between both minima, clearly Ay^^'"-* = 7r/a, while the value of 



^In units of Mp — 1. 



AX^^' ^ cannot be written analytically. If we use this relation between the fields X and 
Y we recover again the standard one-dimensional case where the canonical field will now 

be V^(+'-) = Vi(l + (flSS)') l"^- 

In order for this model to be useful for explaining the actual acceleration stage of 

the universe and also the smallness of the cosmological constant we should check that 

the tunneling probability for the decay from V^''^'* to V^*- ^ is larger than the tunneling 

probability for the decay from Vq^-* to the minimum at F = 0, and also that the minimum 

at Vq ^ has a decay time larger than the life of the universe. 

Let us begin by computing the decay probability from the minimum at yj'*"-' to the 

minimum at F = and comparing it with the decay probability from the minimum at (~ | 

to the minimum at V^'' \ For those cases we will have from that the probability 

will be given in both cases by 

, ^ _ ( 

(1 + 4yJ+V35i2)2 

lat the decay witi 

most probable. The tension of the bubble wall in both cases can be written as 



as V^'''*''* ^ ^- Therefore it is clear that the decay with smaller tension Si will be the 



S- 



.(+,oo) _ ^y{+,oo) ^^(+,oo) ^ < ^^^) A^(+-) . (6.6) 



O 

where Vi'^'°^\ V^'^' denote the height of the maxima that separates any two minima ^ 
(see figure and ^). Also A(^''"*''°°-' , Ay?'^'^'^-' denote the variation of the canonical field i^jj) 
in each case. As can be shown in figure ^ we have that 1 > A99^"^'°°'* ^ A^?'-"*"'"-' and ^ 

we will have in general terms that S'|''''°°'' > S^-^' \ Therefore 
we find that in this case it is more probable to decay to the minimum 

o 

to the minimum V = at X ^ 00. 

Now we must compute the probability of decay from the minimum at ^ to the (J-| 
minimum at V = 0. In this case the tension of the bubble wall can be written as 



5;-'°°)~^f/-'°")Av?(-'°°), (6.7) 

where '"^^ is the height of the maximum that separates the two minima. From the form 
of the potential shown in figure we can assume that Aip ~ 0(1), so it is clear from ( |6.7] ) 
that Si ^ Vq \ Therefore for the decay probability one simply gets 

P(_ ^) ^ exp ^ ~ exp ^ + — ^ + • • • . 6.8 

For V^*- ^ ~ 10^12'^ this probability is extremely small, so therefore, this dS vacuum can be 
considered stable in practical terms. Also note that the thin-wall approximation is always 
true within this model as 0(1) > Aip and p, (V^o''^'*)^^^^ ^ ^• 

It is interesting to note that the analysis of the stability of the last minimum in the 
= 1* is very similar to this last case. The reason for that is the following: the analytical 
discussion of the scalar potential in the non-supersymmetric case developed in section ^ 



-23- 



showed that apart from the minimum we also find a saddle point located in (for D = Dq) 

X° = 1.12\/]V + 1.32 + 1.33^ + ••• (6.9) 



y° = 1.44\/]V + 0.37- 1.36^ + ••• (6.10) 



with a value of the potential given by 

K = 0-11]^ - 0-46^ + 1.57^ + • • • (6.11) 

This value of the potential it is small for large A^, but is still big compared with lO"^^'^ 
for reasonable values of A^. As this point is a maximum in one direction but a minimum 
in the other directions, we can perform an analysis of the stability of the vacua following 
the lines of the two exponential case. As in this case we also find A99 ~ 0(1) and Vsp ^ 
V^jj^ ~ 10"^^*^ we will arrive to the same conclusion as the one in ( |6.8| ). 

In the N = 1* potentials, as in all the cases with many minima, a typical situation 
may be that some of the minima would correspond to de Sitter space and others to anti de 
Sitter. We may imagine living in the one corresponding to de Sitter space with the smallest 
value of the cosmological constant and would wonder about its decay probability towards 



a global minimum with negative cosmological constant. As discussed in |21], the decay to 



a state of negative vacuum energy may or may not occurr and may lead to gravitational 



collapse. This has been recently reanalysed in |24]. 

Finally we would like to mention that it is also necessary to check that the decay 
times of the de Sitter vacua are also not too long. The fact that a de Sitter space has 
finite entropy introduces a time scale that is the Poincare recurrence time tr ||2^. This 
quantity is given by tr ~ e'^'*'^, where Sds denotes the entropy of the de Sitter space. For dS 
space the entropy has a simple sign-reversal relation with respect to the euclidean action 
calculated for the false vacuum dS solution ip = ipQ, which is given by -247rVVb. Then 
the recurrence time can be written as tr ~ e^^"^ Z^". An interesting property of this kind 
of models is that the decay time of the de Sitter vacua never exceeds the recurrence time 
of the de Sitter space tr- This was first noticed in and can be easily checked from the 
following expression 

^^^"^ ln(l + 4yJ±V35?)' >0 =^ tr>tdecay (6.12) 



In ^decay 

The problems related to the decay time fdecay exceeding the recurrence time tr will then 
not appear in these models. 



7. Discussion 



We have presented examples of multiple de Sitter vacua in string theory. Even though 
the examples we consider are still relatively simple, they illustrate what can be expected 
from the general vacuum structure of string theory, i.e. a multitude of vacua with different 
values of the cosmological constant. 



The parameters of the theory allow for one of the minima to have a cosmological 
constant as small as we want. This requires fine tuning^ but it can be ameliorated given 
the large number of minima, indicating an anthropic approach to the cosmological constant 
problem, as advocated by different authors [T^. The number of minima increases 
with the rank of the gauge group N. Furthermore since there is an underlying SL(2, Z) 
symmetry behind the M = 1* theory, we may expect that there could be further, possibly 
infinite, minima if we explore other fundamental domains of this group. We have essentially 
only explored N copies of the strips defined by —1/2 < Y < 1/2 for X outside the unit 
circle in the upper half plane. For each of these strips, the modular group has an infinite 
number of fundamental domains that could indicate a huge multiplication of the number 
of minima, inequivalent from the ones we found since SL(2, Z) is not a symmetry of the 
theory.^ However their study is beyond the limit of validity of our effective actions which 
are trusted only for X greater than the string scale. Furthermore our potentials have 
periodicity in the Y direction. Any correction that would slightly break this periodicity 
could give rise to an infinite number of minima. 

The large number of vacua that can appear in these theories due to the nontrivial 
superpotential complements the already rich structure of vacua due to the presence of the 
fiuxes [^^. The large number of 3-cycles in typical Calabi-Yau manifolds imply a large 
amount of possibilities. Remember also that although the combination of the fiuxes KM 
is restricted by the tadpole cancellation condition, the ration K/M is a free (quantised) 
parameter.^'' This combination appears in the warp factor and allows the tuning of the 
parameter D to get a small cosmological constant, defining the 'discretuum' of vacua as de- 
scribed in 1^, 12, 27 1 . All these effects were present in the single exponential case considered 



in KKLT. The large number of vacua we found has to be multiplied by this degeneracy. Al- 
though it may not be large enough degeneracy for a naturally small cosmological constant, 
the greater the number of minima the more natural is to find a cosmological constant of 
the right size. 

The fact that our potentials depend nontrivially on at least two real fields X, Y makes 
the discussion of the system more interesting than for single field potentials in several 
ways. Since we may have many de Sitter minima, if we imagine the universe starting in 
any of them, it would leave naturally to different periods of infiation, either from tunneling 
between minima but also by naturally rolling after the tunnelings. There are so many 
valleys and hills in the potential that it may not be impossible to find regions of slow roll 
between different minima. 

This combination of tunneling plus rolling has been considered in the past on different 



models of infiation such as open inflation |29, 30 1. A detailed study of the possibility for 



*The fact that we do not need to have a nonvanishing superpotential from the fluxes (Wo) indicates 
we do not need to fine tune this quantity as in [^. However we stiU need to fine tune the supersymmetry 
breaking parameter D which even though is discrete it can be varied almost continuosly j^, 

^Notice also that in computing the scalar potential we restricted ourselves to one single phase of the 
A/" = 1* theory. In general we could consider any other values of p and q. 

^"Notice that typical four-folds can have Euler number between 10'^ and 10® [Q, allowing for many 
different combinations of M and K to satisfy the tadpole condition. 



these potentials to give rise to realistic (eternal) inflation would be clearly of great interest. 
For this we recall that the main obstacle to have successfuU inflation from string theory is 
precisely the lack of control of the moduli potentials. In particular the different proposals 
of D-brane inflation |31] assume that there is an unknown stringy mechanism that fixes the 
moduli and then, after that, D-brane inflation could occur. This is clearly a very strong 
assumption and so far attempts of combining D-brane inflation with moduli fixing have 
been running into problems. Therefore we may consider seriously the possibility that 
actually the modulus T could be the infiaton field, and these potentials could give rise to 
interesting combinations of inflationary processes. This is clearly a possible subject for 
future investigation. 

Another open question left unanswered here is the detailed analysis of the decays of 
the different minima. This is complicated by the facts that we have two-field potentials 
and that the fields do not have canonically normalised kinetic terms. Our discussion 
was mostly carried for the simplest case of two exponentials with only two minima, but 
clearly a complete analysis for the M = 1* and other more general potentials would be 
needed. 

Even though our models are relatively simple, they illustrate the potential richness of 
the string landscape once the different 'moduli' fields acquire nontrivial potentials. We do 
not pretend that the superpotentials discussed here, such as the one from the N = 1* the- 
ory, would be particularly special over other more realistic realisations of non-perturbative 
potentials. Actually, we regard it as an interesting tool which includes all the ingredi- 
ents expected from non-perturbative physics, in particular the infinite instanton sum can 
be under control thanks to the mathematical properties of modular forms. Other non- 
perturbative superpotentials recently derived, including different deformations of = 4 
theory [^] would be interesting to study. 

The mechanism of KKLT, although very interesting, takes the simplest class of models 
in which only one modulus is left unfixed by the fiuxes. In more realistic models we would 
expect many moduli left unfixed by the fiuxes and finding a many fields non-perturbative 
potential for them would be an important challenge. We regard our study here as a 
nontrivial, yet manageable, realisation of the multi-vacua scenario in string theory and 
hope it may be of interest to explore further properties of the space of string vacua. 



8. Note added in revised version 

The analysis presented in this article was performed using a supersymmetry breaking term 
induced in the potential due to the inclusion of an anti D3-brane in the configuration. Such 
a term was taken to be given by 5V = D / , with D a positive constant depending on the 
fiuxes turned on in the compactification. Nevertheless, when this work was already finished 
it was pointed out in [^2| that this induced term is not the one to be used if the anti D3 brane 
is in the warped region, as it is. The correct term should be 5V = D/X"^ (with D also a 
positive, fiux-dependent, constant). This arises because the supersymmetry breaking term 

^^S. Kachru and J. Maldacena private communication ^A- 
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Table 7: Minima of the scalar potential for several non-supersynimetric cases. 



coming from the inclusion of an anti D3-brane in the warped compactifications considered 
here scales like e^'^/X^, and in the highly warped regime e^^ ~ Xe"^'^^^'^^"^'^ . This fact 
does not change the main conclusions of the paper although it changes the numerical 
analysis. We briefly summarize in this added note the changes produced in the numerical 
analysis. 

As it was found in the D/X^ case, with the introduction of the supersymmetry breaking 
term D / X'^ it is also possible to lift the vacua from anti de Sitter vacua to de Sitter vacua, 
for some range of the parameter D. Again we found that the effect of the supersymmetry 
breaking term depends on the range of values of the parameter D. If D is very small 
(compared with the value of the potential in the supersymmetric case), the potential will 
not change substantially and the minima remain anti de Sitter. For a critical value of D one 
of the minima will move up to zero vacuum energy and then to de Sitter space. Continuing 
increasing D, more minima become de Sitter until all of them are either de Sitter or 
Minkowski. We will denote by Dq that precise value of the parameter D such that all the 
values of the potential at the minima are positive or zero. If D > Dq the nonsupersymmetric 
term starts dominating the potential and starts eliminating the different extrema to make 
the potential runaway with X. 

Actually, we have computed numerically the minima of the potential for D = Dq for 
several values of A'^, such as = 50, = 100 and A'^ = 500. The results obtained from the 
analysis are shown in table ^, where in those cases Dq = 0.04836, 0.04958, and 0.04969 for 
N = 50, 100, and 500 respectively. 



Also we show in figure |9| the scalar 
potential for N = 100 when the super- 
symmetry is broken. Note that com- 
paring with figure |^ the main differ- 
ences are that now the potential is less 
smooth and therefore has more minima 
(this can also be seen by comparing the 
information shown in table |5| and ta- 
ble 0). 

In table ^ we show a similar in- 
formation as the one shown in table ^ 
where we write the number of minima 
varying with N. In table ^ we denote by 
Minx,y the minima with largest value 
of X, Y and by Minv'=o the minima 
with vanishing value of the potential. 

We can notice that the number of minima increases with N more or less in the same rate 
than in the supersymmetric case and then faster than in the D/X^ case. The reason for 




Figure 9: Graph for the non-supersymmetric scalar 
potential with N = 100. 
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Table 8: Minima for different values of N in a non-supersymmetric case with D = Dq. 
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Figure 10: Potential versus X ioY N — 100 and different values of D. 



this is that the value of D for which all the anti de Sitter minima of the potential get lifted 
is smaller than in the D/X^ case, and therefore in the present case the 5V term smooths 
less the potential. 



Also, in figure IC we illustrate the effect of the corrected non-supersymmetric term in 
the potential. The value of the potential at the minima is presented for several values of the 
parameter D. For = we have the super symmetric case, the increasing of the parameter 
D will reduce the number of minima and increase the value of the compactification scale 
at the minima. Note that now the minima are not ordered for increasing X, and also that 
the values of the compactification scale at the minima are, in general, smaller than in the 
D / case but still large enough so that our approximations remain valid. 

Summarising: the difference with the 1/X^ potentials is that now there are more de 
Sitter minima, they are no longer ordered with respect to the size of the extra dimension 
and the value of the volume at the minima tends to be smaller than in the 1/X^ case, but 
still large enough for the large radius approximation to be trustable. 

Finally we would like to point out that the 1 / X^ behaviour of the nonsupersymmetric 
part of the potential presented in the main text is still relevant. This is due to the recent 
proposal in [5^] in which instead of introducing an anti D3 brane on the throat, a flux of 



magnetic fields on the D7 branes is considered. This adds an extra term in the potential 
which happens to be proportional to if the D7 brane is not sitting on the throat 

and to l/X"^ if the D7 brane is on the throat. In each case the rest of the analysis is as 
presented in here and therefore both potentials are relevant depending on the location of 
the D7 branes. 
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